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ABSTRACT 

We report on the detection of giant pulses from the Crab Nebula pulsar at a frequency of 200 MHz using 
the field deployment system designed for the Mileura Widefield Array's Low Frequency Demonstrator 
(MWA-LFD). Our observations are among the first high-quality detections at such low frequencies. The 
measured pulse shapes are deconvolved for interstellar pulse broadening, yielding a pulse-broadening 
time of 670±100 (is, and the implied strength of scattering (scattering measure) is the lowest that is 
estimated towards the Crab nebula from observations made so far. The sensitivity of the system is 
largely dictated by the sky background, and our simple equipment is capable of detecting pulses that are 
brighter than ~9 kJy in amplitude. The brightest giant pulse detected in our data has a peak amplitude 
of ^50 kJy, and the implied brightness temperature is 10 316 K. We discuss the giant pulse detection 
prospects with the full MWA-LFD system. With a sensitivity over two orders of magnitude larger than 
the prototype equipment, the full system will be capable of detecting such bright giant pulses out to a 
wide range of Galactic distances; from ~8 to ~30 kpc depending on the frequency. The MWA-LFD will 
thus be a highly promising instrument for the studies of giant pulses and other fast radio transients at 
low frequencies. 

Subject headings: pulsars: general - pulsars: individual (Crab pulsar) - ISM: structure - ISM: 
individual (Crab Nebula) - scattering 



1. INTRODUCTION 

Giant pulses from the Crab Nebula pulsar have been 
known since the early days of pulsar observations (Staelin 
& Reifenstein 1968; Argyle & Gower 1972). These spo- 
radic, large-amplitude bursts c an exceed mean pulse ener- 
gies by factors as large as 10 4 (jLundgren et al.lll995[ ), and 
they remain as an intriguing phenomenon for theories of 
the pulsar emission mechanism to explain. Their short du- 
rations also allow them to serve as sensitive probes of the 
nebula and the interstellar medium. Much has been learnt 
from observations so far, and it is fairly well established 
that: (i) the fluctuations in their amplitudes are associ- 
ated with changes in t he coherence of the radio emission 
(jLundgren et al.l [l995). (ii) the emission i s broadband - 
extending over several h undreds of MHz (jSallmen et al.l 
119991: IPopqv et all I2006D . and (iii) they are superposi- 
tions of extremely na rrow nanosecond-duration structures 
(Hankins et al. 20Q3f). The Crab giant pulses have some of 
the largest implied brightness temperatures of known as- 
trophysical sources, and are known to reach 10 32 K at mi- 
crosecond time resolutions (jCordes et al.l 12004), and even 
highe r (10 37 K) at nanosecond resolutions (|Hankins et ail 
120031) . 

Although the Crab pulsar was o riginally discovered at 
a low radio frequency of 111 MHz (Staelin & Reifenstein 



Il968h through detection of its giant pulses, the bulk 
of the observational studies of it over the past three 
decades has focuse d primarily at higher frequ encies (e.g. 
iMoffett fc Hankinsl 119961 : ICordes et all 120041 ). This is 
partly due to limitations imposed by the dispersive and 
scattering effects in the interstellar medium (ISM) that 
degrade the signal strength at low frequencies. A more 
significant reason perhaps is that most existing large in- 
struments are optimised for best performance above 300 
MHz. For single dishes this has been driven by a quest 
for sensitivity, and for interferometers by resolution and 
sensitivity, and also the need to avoid the problems due 
to ionospheric distortion of astronomical signal at low fre- 
quencies. The recent enormous increase in affordability of 
computing and the impressive advances in the capabilities 
of digital hardware have brought the computational needs 
of low frequency radio imaging within reach. This has led 
to a resurgence in the interest in low frequency astron- 
omy and a number of instruments like the Mileura Wide- 
field Array Low Frequency Demonstrator (MWA-LFD) in 
Western Australia, Low Frequency Array (LOFAR) in the 
Netherlands and Long Wavelength Array (LWA) in South- 
west US are currently in varying stages of design, devel- 
opment and deployment. The ever growing problem of 
man made radio frequency interference (RFI) continues 
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to remain a major concern however, especially in densely 
populated parts of the world. 

Despite the paucity of low- frequency studies of the Crab, 
a wealth of information can be gleaned from such studies. 
As such pulsar radio emission is strongly frequency depen- 
dent, with most intrinsic and extrinsic effects becoming 
more pronounced at lower radio frequencies. In partic- 
ular, the Crab pulsar is well known for its highly pecu- 
liar and complex fr equency evolution of its pulse morphol- 
ogy a nd structure ( Cordes et alJl2004l iMoffett k, Hankinsl 
1996). Therefore it would be interesting to examine how 
the properties of pulse emission evolve at low radio fre- 
quencies. Further, the scattering and dispersion effects 
also become more prominent with decreasing observing 
frequency, and are best studied at low frequencies. 

The MWA-LFD is a promising instrument with capabili- 
ties that will complement those of existing large telescopes, 
particularly the frequency coverage. Given its high sen- 
sitivity, multi-beaming and wide field-of-view capabilities, 
this instrument promises excellent prospects for pulsar and 
transient science. In this paper, we report on successful 
detection of the Crab giant pulses at 200 MHz using the 
prototype system that was built for its early deployment 
activities. Although the Crab pulsar was detected earlier 
at a nearby frequenc y (196.5 MHz) via its average emission 
(jRankin et al.lll970l ). our observations mark the first de- 
tection of individual giant pulses at 200 MHz, and the data 
are used to estimate the scattering due to the intervening 
ISM as well as the pulse amplitudes and brightnesses. We 
also assess the detectability of giant pulses with the full 
MWA-LFD system. 

2. OBSERVING SYSTEM AND DATA PROCESSING 

2.1. The MWA-LFD Early Deployment system 

The MWA-LFD array design covers a broad frequency 
range from 80 to 300 MHz, and features a large number 
of small phased-array antenna systems or tiles, each con- 
sisting of a 4 x 4 array of dual-polarisation wideband ac- 
tive dipoles and an analog delay- line beam- former unit. 
These tiles yield a field of view of — 200A 2 deg 2 (where A 
is the observing wavelength in meters), i.e., approximately 
20° in angular extent at a frequency of 200 MHz. Three 
such tiles were deployed at the Mileura station in Western 
Australia (26 26'S, \VJ°12'E) as part of the "early deploy- 
ment" (hereafter ED) program, with baselines of «146 m, 
«325 m and «212 m. These tiles, plus a simple data cap- 
ture and software correlation system, form the observing 
system used for our data collection. Further details on the 
ED system and activities are described in Bowman et al. 
(2006). For an overview of the array design and science 
goals see http://haystack.mit.edu/ast/arrays/mwa. The 
full system, when operational, will consist of 512 such an- 
tenna tiles, more than 95% of which will be distributed in 
an area — 1.5 km in diameter. The remaining tiles will be 
distributed in a concentric ring with an outer diameter of 
-3.0 km. The entire band from 80 to 300 MHz will be 
sampled, 32 MHz of which will be processed and recorded. 

2.2. Observations of the Crab Pulsar 

Observations of the Crab pulsar were made on 22 
September 2005, during the last of four two- week long 
expeditions conducted as part of the ED program. Data 



were recorded over a bandwidth of 8 MHz (approximately 
6 MHz of useful bandpass) centred at frequencies 100 and 
200 MHz for a total duration of approximately 3.5 hours 
each. The signal was dig itally sampled at 1 6 MHz using 
the "Stromlo Streamer" (jBriggs et al.ll2007ft using 8 bits 
and then re-encoded into 4 bits using a scheme similar to 
"//-law" companding. This scheme uses increasingly larger 
spaces between quantisation levels which allows the input 
power to vary by several orders of magnitude without 
incurring significant losses in sensitivity. 

Four independent data streams were recorded at a given 
time. The system was configured to record one of the two 
orthogonal polarizations from two tiles, and both polar- 
izations from the third tile. This implied an aggregate 
data rate of 32 mega bytes per second, i.e., a raw data 
volume of «0.8 tera bytes. The data were stored on disks 
for offline processing and all processing operations, such as 
detection and de-dispersion, as well as generation of spec- 
trometer (filtcrbank) data and phased-array data streams, 
were performed in software using the Swinburne super- 
computing facility. 

2.3. Data Processing and Analysis 

Data processing was performed in multiple stages. The 
initial processing was done using software that was devel- 
oped for calibration of the ED system. This was tailored 
to output an incoherent filterbank stream of data at a 
time resolution of 1024 (is and a spectral channel reso- 
lution of 8 kHz, extracting only the central 4 MHz part 
of the band. This coarse resolution processing, which re- 
sults in significant dispersive smearing («471 /xs), served 
as the preliminary processing pass in order to examine it 
for visible signatures of giant pulses before embarking on a 
more detailed and rigorous analysis. Though it limited our 
sensitivity to only the brighter giant pulses, it still proved 
very useful. 

The data were then subjected to a blind search for giant 
pulses by performing standard procedures such as dedis- 
persion and progressive smoothing of the dedispersed time 
series, and identifying samples that are above a set thresh- 
old. Dedispersion was performed at multiple trial values 
around the Crab's nominal dispersion measure of 56.791 
pc cm -3 in order to confirm the detections. This analysis 
yielded a handful of confirmed detections at 200 MHz and 
none at 100 MHz. The data at 200 MHz were subsequently 
reprocessed to output a filterbank stream at a much higher 
time resolution of 256 /is and an improved spectral reso- 
lution of 4 kHz. A similar giant pulse search procedure 
was then repeated. The implied dispersive smearing in 
this case (235 fjs) is quite comparable to the time reso- 
lution. This procedure confirmed all the detections from 
the first coarse-resolution search, but did not result in any 
additional detections. 

In order to overcome the inherent limitations in process- 
ing such incoherent filtcrbank data, we developed software 
that directly operates on raw voltage samples and per- 
forms a phased-array addition of the signals, eventually 
constructing a coherent filterbank stream of data0. Details 
i 

The dispersive delay, t(j m f» 8.3 fis DM Au/u 3 where DM is the 
dispersion measure, Av is the spectral channel width (MHz) and v 
is the frequency (GHz) of observation. 
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of this scheme are described in the following two sections. 

2.3.1. Generation of Phased- array Data 

The underlying principle here is similar to that of the 
tied array mode of operation as employed in instruments 
such as the Australia Telescope Compact Array (ATCA) 
and the Westerbork Synthesis Radio Telescope (WSRT), 
or the phased-array mode at the Giant Metrewave Radio 
Telescope (GMRT). Here, voltage samples from individual 
array elements are summed after correcting for geometric 
and instrumental delays in arrival times of the signals. For 
our ED system, the combination of rather short baselines 
and long wavelengths means such corrections are largely 
dominated by the array layout or geometry, rather than 
any tile-dependent phase corrections as derived from the 
calibration process. In fact, applying the latter results 
in only a marginal improvement in the signal strength. 
While instruments such as the ATCA and the WSRT are 
equipped with dedicated hardware to perform such an op- 
eration, in our case, this phased-array addition of signals 
was fully realized in software. 

For our 3-tile ED system, we would expect the signal to 
be improved by a factor of y3 in this manner, but in prac- 
tice even more by performing a coherent dedispersion of 
the signal. Ideally, the phased array summation will need 
to be performed separately on each of the two orthogonal 
polarizations, however in our case only the east-west polar- 
ization has signals recorded for all three tiles. The phased- 
array data stream after the summation is simply written 
out as floating point numbers, mainly for the simplicity 
of the data format and compatibility with the processing 
software. This resulted in a factor of four increase in the 
net data volume («1.5 tera bytes per frequency). The data 
were transfered onto disk drives and made accessible to the 
Swinburne supercomputer via the Grange-Net facility. 

2.3.2. Dedispersion and Giant Pulse Search 

Our phased-array data streams were written out in the 
form of voltage samples (i.e., baseband data format), and 
hence allowed a coherent dedispersion to be performed on 
the signal. This procedure removes the dispersion entirely, 
and is especially important at low radio frequencies where 
the signal degradation due to dispersion can be rather 
severe. The dedispersed data are then subjected to the 
pulse detection algorithm to search for giant pulses. The 
procedure s we have adopted a re very similar to that de- 
scribed in iKnight et al.l |2006) . The voltage samples are 
first Fourier transformed to the frequency domain and the 
spectra are divided into a series of sub-bands. Each sub- 
band is then multiplie d by an inverse response filter (ker- 
nel) for the ISM (e.g. lHankins k Rickettjll975r ). and then 
Fourier transformed back to the time domain to construct 
a time series with a time resolution coarser than the orig- 
inal data. By splitting the input signal into several sub- 
bands, the dispersive smearing is essentially reduced to 
that of an individual sub-band. This also means the proce- 
dure uses shorter transforms than single-channel coherent 
dedispersion and is consequently a more computationally 
efficient method. 

This "coherent filterbank" stream of data are then 
square-law detected and combined after correcting for dis- 

Data are coherently dedispersed prior to forming a filterbank stream. 
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Fig. 1. — A bright giant pulse detected with the MWA-LFD early 
deployment system. Top: Dedispersed time series of a short data 
segment near the giant pulse. Middle panels: Diagnostic plots from 
the giant pulse search processing; plots of signal-to-noise ratio vs 
dispersion measure and time resolution, and dedispersed time series 
for 8 sub-bands, each 1 MHz wide, within the 8 MHz bandwidth. 
Bottom: The dedispersed pulse shown in high time resolution. The 
scattering tail is clearly visible. 



persive delays between the sub-bands, to construct a single 
coherently dedispersed time series for the entire band. The 
pulse detection procedure involved progressive smoothing 
(convolution) of time series with matched filters of widths 
ranging from 64 /is to 2048 /is in steps of a factor two, and 
identifying the intensity samples that exceed a set thresh- 
old (e.g. 8.5 a). In addition to performing dedispersion 
at the Crab's nominal DM of 56.791 pc cm" 3 , we also 
perform this procedure over a large number of adjacent 
DM values (typically over a DM range «0.1DM cra b, in 
steps of 0.01 pc cm -3 ), in order to confirm the dispersion 
signature. For each DM and the matched filter width, we 
compute signal-to-noise ratio (S/N) of the pulse amplitude 
over a short stretch of time centered on the pulse. This 
process is computationally intensive and was carried out 
on Swinburne's supercomputer. From this analysis, diag- 
nostic plots are generated for each candidate giant pulse 
as shown in Fig. [TJ We also check if the signal is broad- 
band by displaying the dedispersed time series in 8 x 1 
MHz sub-bands. These diagnostic plots were subjected to 
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a careful human scrutiny to discriminate real giant pulses 
from spurious signals. 

2.3.3. System Sensitivity and Flux Calibration 

The Crab Nebula is a fairly bright and extended source 
in the radio sky, with a flux density of ~955 is~ - 27 Jy 
(|Bietenholz et al.1ll997l ) and a characteristic diameter of 
~5'.5. Thus, in general, for most conventional observa- 
tions with large single-dish telescopes such as Arecibo and 
Parkes, the nebul ar emission dominat es the system noise 
temperature (e.g. ICordes et al.1l2004h . However, in our 
case the combination of a low antenna gain, large sky 
background, high system temperature, large field of view, 
and significant side lobes of the antenna response function, 
leads to a very different scenario whereby the system tem- 
perature T sys is heavily dominated by the sky background 
itself. This is further compounded by the fact that the tile 
response pattern becomes more complex at large zenith 
angles that are relevant for the Crab, and varies signifi- 
cantly over the duration of the observation depending on 
the zenith angle (ZA). Further technical aspects of the ED 
system are discussed by Bowman et al. (2006). 

In order to account for these effects, and to obtain a re- 
alistic estimate of the sky background temperature, T s k y , 
we convolve the sky model, T s k y (6), with the antenna 
power response function, W(6>). The sky model T s k y (9) 
is obtained by extra polation of the 408 MHz sky map of 
lHaslam et al.l ( 19821 ) assuming a spectral index a — —2.6. 
For our purpose, W(0) is the phased-array beam pattern 
when the tile beam is pointed toward the Crab. Perform- 
ing the above calculation yields T s k y = 180 K when the 
Crab is at transit (ZA = 48°), but may increase as much as 
by 20% at larger zenith angles of the Crab (e.g. ZA = 60°). 
The net system temperature is given by T sys = T s k y + 
T reC v , where T recv is the receiver temperature. For the ED 
system, the receiver temperature is a strong function of fre- 
quen cy, and is 200 ± 20 K near 200 MHz (jBowman et al.l 
120061) . Thus T sys = 380±30 K at 200 MHz, taking into 
account the nominal uncertainty in receiver temperature 
measurements and change in T s k y with the ZA. 

In order to translate our system temperature estimates 
to their equivalent flux densities (i.e., the system equiva- 
lent flux density, S sys ), we require an estimate for the array 
gain, G, expressed in conventional units of K Jy -1 . This 
gain parameter can be derived from the effective area of 
the tiles, A e fi , using the simple relation, G = A c g- /2 ke , 
where ke is the Boltzmann constant. For the MWA-LFD 
antenna, the effective area is a strong function of frequency 
and is given by A c g = A 2 /x, where x is the frequency- 
dependent scaling parameter. Based on our design simu- 
lations x = 2.1 at 200 MHz, i.e., A off w 1 m 2 for a single 
dipole. Thus for the ED system which comprises 48 dipoles 
in total, the net effective area is w 50 m 2 at 200 MHz, i.e., 
an equivalent system gain of «0.02 K Jy -1 . However, as 
described in Bowman et al. (2006), the tile gain also has a 
strong functional dependence with the zenith angle. The 
measured dependence of cos 2 6 (Z A) is quite steeper than a 
theoretically expected cos 2 (ZA) form, and thus the nomi- 
nal gain for the ED system is «0.008 KJy -1 towards the 
Crab. Further, the zenith angle varies from 60° (near rise 
or set) to 45° (at transit) for the Crab, leading to a gain 
variation by a factor 2.5 over the observation. Thus, our 
T sys estimate of 380 K corresponds to an equivalent sys- 
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Fig. 2. — Measured and reconstructed pulse shapes of a bright 
giant pulse at 200 MHz. The best fit scattering function (i.e., the 
pulse broadening function, PBF) is shown in the inset. The top 
and bottom panels correspond to different models of the PBF; an 
exponential PBF that corresponds to an infinite thin screen model 
of scattering medium (top), and a more rounded PBF that cor- 
responds to a uniform distribution of scattering medium (bottom). 
The deconvolution results in a significant increase in the peak signal 
strength and a pulse shape that is presumably closer to intrinsic. 



tern noise of S sys = T sys /G » 47.5 kJy. For our recording 
bandwidth of 8 MHz and assuming a typical matched filter 
width of 250 [is, the effective S sys is of the order of ~1100 
Jy, i.e., a minimum detectable pulse amplitude of ~5.5 kJy 
for a 5-ct detection threshold. 

2.3.4. Deconvolution of Pulse Broadening 

Fig. [T] shows an example of giant pulse detection. As 
evident from this figure, the measured pulse shape is quite 
asymmetric and is marked by a clear lengthening at the 
trailing end, a characteristic of the pulse broadening phe- 
nomenon that results from multipath scatt ering in the 
intervening interstellar med ium (ISM) (e.g. IWilliamsonl 
Il972t ICordes fc Laziol |2001|) . This effect is more pro- 
nounced in obse rvations of dista nt pulsars and at low 
frequencies (e.g. iBhat et al.l l2004f ). The measured pulse 
shape, y(t), can be modelled as effectively a convolution of 
the intrinsic pulse shape, x{t), with the impulse response 
of the scattering medium, g(t), often referred to as the 
pulse broadening function (PBF). More precisely, 

y(t) = x{t) ® g{t) <g> r(t) (1) 

where r(t) is the net instrumental resolution function. The 
observed pulse broadening is usually quantified by a time 
scale, Td , characteristic of a PBF fit to the measured pulse 
shape. The exact form of the PBF and its scaling with 
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frequency depend on the nature of distribution of scatter- 
ing material al ong the line-of-sight and o n its wave number 
spectrum (e.g. lLambert fc Ricketdll999l ). 

In order to estimate the pulse broadening, we adopt 
the CLEAN-based deconvolution approach developed by 
iBhat et"ail (|2003f ). Unlike the traditiona l freq uency- 
extrapolation approach (e.g. lLohmer et alj l200lh . this 
method makes no prior assumption of the intrinsic pulse 
shape, and thus offers a more robust means of determin- 
ing the underlying PBF, and therefore can be applied to 
a wide variety of pulse shapes and degrees of scattering. 
The procedure involves deconvolving the pulse shape in 
a manner quite similar to the CLEAN algorithm used in 
synthesis imaging, while searching for the best-fit PBF 
and recovering the intrinsic pulse structure. It relies on 
a set of figures of merit that are defined in terms of pos- 
itivity and symmetry of the resultant deconvolved pulse 
and some parameters characterizing the noise statistics in 
order to determine the best-fit PBF. Such an approach is 
especially justified for the Crab giant pulses, as they are 
known to show structures at timescales down to microsec- 
onds or even nanoseconds, and show complex and quite 
unusual evolution in pulse morphology with frequency. 

The deconvolution procedure e mploy ed here is similar 
to that described in IBhat et all (|2003l ) , except that the 
restoring function used in our case is a much simpler one. 
This is because the pulse smearing due to residual dis- 
persion and instrumental effects is negligible. The major 
contributing factor to our restoring function is that due to 
time binning of the pulse profile. We used two different 
trial PBFs; the first, gi , is appropriate for a thin slab scat- 
tering screen of infinite transverse extent (PBFi ), and the 
second, g2 , corresponds to a uniformly distributed scat- 
tering medium (PBF2). Their functional forms are given 
by (e.g. lWilliamsonlll972D . 

g 1 {t) = T- 1 cxp{-t/T d )u{t), (2) 

g 2 (t) = (TrVj/^ 5 ) 1 / 2 exp(-vr 2 T d /At) u(t). (3) 

where u(t) is the unit step function, u(< 0) = 0, u(> 0) = 
1. While PBFi , which has a one-sided exponential form, 
has been commonly used by several authors in the past, 
the latter PBF is a generic proxy for more realistic dis- 
tributions of scattering material. The best-fit PBFs ob- 
tained in this manner along with the respective intrinsic 
pulse shapes are shown in Fig. [5] As evident from this 
figure, the deconvolution process naturally leads to a sig- 
nificant increase (nearly by a factor two) in the effective 
signal-to-noise ratio of the pulse. 

3. RESULTS AND DISCUSSION 

3.1. Giant Pulse Amplitudes and Brightness 
Temperatures 

The measured amplitudes of giant pulses at low radio 
frequencies are determined by a number of factors such 
as the nature of the emission spectrum, the degree of 
scattering and the time resolution of the data. At fre- 
quencies £ 300 MHz, the pulse energies are known to 
follow a i^~ 4 - 5 form, i.e., a spectral form that is much 
steeper than typical pulsar spectrum. The behaviour of 
the spectrum at low frequencies is not well understood, 
although some of the early work suggests a turnover of 
the spectrum near our observing frequency of 200 MHz 
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Fig. 3. — Summary plot of giant pulse detections at 200 MHz. 
The peak flux densities range from 10 to 50 kjy. Despite small- 
number statistics, these detections seem to follow a roughly power- 
law distribution in pulse amplitude and energy. In total, 31 pulses 
were detected in 3.5 hours of observations. The S/Ns shown are 
for the east-west polarisation; the signal is detected in north-south 
polarisation in most cases, albeit with a poorer S/N. 



(e.g. iManchester fc Tavlorl Il977r j . Further, although the 
Crab pulsar is known to show quite a peculiar evolution in 
frequency in terms of pulse morphology and structure of 
both the average and giant pulse emission, at frequencies 
& 300 MHz, scattering is likely to dictate the measured 
pulse shape as the pulse smearing becomes severe. Thus, 
even though the giant pulse emission inherently possesses 
structures down to microsecond or even nanosecond levels, 
finer structural details of individual giant pulses are simply 
not discernable due to severe broadening of the intrinsic 
pulse expected at these frequencies. Although in principle 
our deconvolution approach is useful for potentially ex- 
tracting the information in such situations, unfortunately 
the signal-to-noise of our data does not seem large enough 
to attempt such an analysis. The effective time resolu- 
tion of the data is also an important factor, as any level 
of smoothing (convolution) in time performed to aid ei- 
ther a detection of pulse, or an enhancement in the signal 
strength, would mean an inevitable reduction in the peak 
pulse amplitudes. Thus, in short, measurements of giant 
pulse amplitudes and fluxes at low frequencies should be 
treated with some caution, especially when compared with 
similar measurements made at higher frequencies. 

Results from our giant pulse search are summarised in 
Fig. [3] where the measured pulse amplitudes are plotted 
against their times of occurrence. The estimates of S/N are 
converted to units of Jy using the flux calibration proce- 
dures detailed in § 12.3.31 Note that these are the measured 
peak pulse amplitudes, i.e., no correction is applied for the 
enhancement in signal strength achievable by the deconvo- 
lution of scattering. Despite the small-number statistics, 
the pulse energies seem to follow a roughly power-law type 
distribution (index = — 1.7±0.1). For the nominal system 
temperature and gain estimates as described in § 12.3.31 
(380 K and 0.01 K Jy -1 respectively), and given our typ- 
ical processing parameters (a net bandwidth of 6 MHz 
and a nominal pulse width of 300 /is), the effective system 
noise turns out to be of the order ~1100 Jy. This means 
our system is sensitive to pulses brighter than ~9 kjy, as- 
suming a 8.5-cx detection threshold. Our search analysis 
has yielded a total of 31 giant pulses above this threshold. 
The brightest giant pulse detected in our data has a peak 
flux density of 50 kjy and and a width «300 /is. The mea- 
sured amplitudes can be translated to implied brightness 



6 



Bhat et al. 



temperatures (Tb ) using a simple relation (based on the 
light-travel size and ignoring relativistic dilation), 



2k T 



vAt 



where S„ is the peak flux density and At is the pulse width; 
ke is the Boltzmann constant and D corresponds to the 
Earth-pulsar distance. Assuming a peak pulse ampli- 
tude of 100 kJy (i.e., the amplitude of the reconstructed 
pulse in Fig. [2} a nd a nominal distance of 2 kpc to the 
Crab Nebula (e.g. ICordes fc Laziol l2002h . we estimate a 
brightness temperature of 10 316 K for our brightest gi- 
ant pulse. However this is likely to be an underestimate 
of the true brightness temperature given the degradation 
of the signal strength expected due to pulse broadening. 
Interestingly, this estimate is comparable to that derived 
from the bri ghtest giant pulse d etected at 430 MHz us- 
ing Arecibo ( Cordes et al.l 120041 ) where the scattering is 
at least an order of magnitude smaller. Indeed much 
higher brightness temperatures (10 37 K) have been re- 
ported from observations at 5.5 GHz (where the scattering 
is negligibly small) m ade at nanosecond time resolutions 
(jHankins et alJl2003h . Despite the various factors that de- 
grade the signal strength at low frequencies, our observa- 
tions confirm fairly high brightness temperatures for the 
giant pulses at 200 MHz. 

The amplitude of the brightest giant pulse detected in 
our data is compared with similar estimates from recent 
observations in Fig. 01 The plot spans a wide range of fre- 
quencies from ^20 MHz to ~9 GHz, i.e., nearly 3 orders 
in magnitude. Observations at most of the higher frequen- 
cies are from an intensive multifrequency campaign with 
Arecibo (jCordes et al.ll2004[ ). while t hose at 23, 111 an d 
600 MHz are from the recent work of iPopov et al.l ( 20061 ) . 
From our non-detection at 100 MHz, we derive an upper 
limit of ~5 kJy for the bright giant pulses detectable at 
this frequency, by assuming a simple v scaling of our 
measured T s k y and nominal values for the anticipated sys- 
tem sensitivity. This is much larger than the 0.7 kJy value 
for th e brightest observed pulse reported by IPopov et al.l 
(2006). However, this discrepancy can be attributed to the 
fact that it is based on a much shorter observing duration 
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Fig. 4. — Peak flux densities of the brightest giant pulses (in 
1 hr) over a wide range of frequencies from 23 MHz to 8.8 GHz 
(solid symbols). The symbols indicat e data from differen t sources - 
circles denote Arecibo observations of Cordes et al. (2004), triangles 
are recent low frequency observations of IPopov et al.l (120061') . and 
squares denote the present work. The unfilled symbols near 100 
MHz indicate limits derived from a non-detection with MWA-LFD 
(square) and the brightest giant pulse in 15 minutes of observation 
(triangle) . 



(15 min) and hence may not have been quite sensitive to 
the tail end of the giant pulse power-law distribution. It is 
still interesting that the measured peak amplitudes of the 
brightest giant pulses in one hour follow roughly power-law 
trends with frequency, with an empirical relation of the 



forr 



-4.3±0.2 



at frequencies £ 300 MHz (i.e., comparable 



to the nominal spectral slope of giant pulse energies), and 
iy +2.7±o.i a j. i ower frequencies (if we exclude the 111 MHz 
value) . While it is harder to extract any meaningful infor- 
mation about the nature of the spectrum from these mea- 
surements alone, this plot may still serve as a useful guide 
to assess the detection prospects of bright giant pulses. 
Characterization of the spectrum requires reliable mea- 
surements of average giant pulse fluxes or energies, which 
are not easily obtainable for all these frequencies. Clearly, 
scattering will strongly influence the giant pulse detection 
prospects within the MWA-LFD's frequency range. We 
discuss the detection prospects for the full system in § |U 

3.2. Implications of Measured Pulse Broadening 

Our pulse deconvolution yields a Td of 670 ± 100 /xs for 
th ePBFgi , and 395 ±50 /js for the PBF g 2 . As described 
in iBhat et alJ (|2004h . in principle, the figures of merit of 
the CLEAN deconvolution technique can be used to assess 
which PBF better fits the data. Specifically, their param- 
eter, f r , which is a measure of positivity, can be treated 
as a useful indicator of "goodness" of the CLEAN sub- 
traction. We expect / r S 1 for successfully deconvolved 
pulses, while larger values imply slightly overCLEANed 
pulses. For the results in Fig. [21 we obtain values of 1.3 and 
0.85 respectively for deconvolution with PBFi and PBF2 , 
which would mean PBF2 is probably a better fit to the 
data. As PBF 2 represents a scattering geometry with dis- 
tributed material, this suggests that much of the observed 
scattering is possibly due to the intervening ionized gas 
along the Earth-pulsar line-of-sight. 

Our Td estimates for the two PBFs differ roughly by a 
factor two, and this is because Td has different meanings 
for the two scenarios. For PBFi , Td is both e^ 1 point of 
the distribution and the expectation value of t, whereas 
for PBF2 , Td is close to the maximum of the distribution, 
which is at (7r 2 /10)Td = 0.99 Td , while the expectation 
value of t is (7r 2 /2)Td = 4.93 Td (see eqn. [3]). 

The measured pulse broadening can be used to infer the 
scattering measure (SM), defined as SM = dsC 2 (s), 
where C, 2 is the spectral coefficient of the wavenumber 
spectrum of electron density irregularities in the ISM 
and D is the Earth-pulsar distance. It signifies the to- 
tal amount of scattering along the line-of-sight and de- 
pends critically on the nature of distribution of the scat- 
tering material. It can be related to pulse broadening 
via the relation T d « 1.1 W T SM 6/5 z/- 4 - 4 L>, where v is 
in GHz, D is in kpc, and W T is a geometric factor that 
depen ds on the line-of-sight di stribution of scattering ma- 
terial (ICordes fc Rickettl ll998). Using this relation, and 



assuming W T — 1 , we estimate the effective SM for a uni- 
form medium as (1±0.3) x 10~ 3 kpc m~ 20 / 3 . This can be 
compared to a value 1 .7 x 10~ 3 kpc m~ 20 / 3 as predicted 
by the NE2001 model ([Cordes fc Lazidl2002ft . 

While the estimates of pulse broadening time are, in 
general, reliable measures of the integrated scattering 
due to the ISM, some caution is needed while interpret- 
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ing the measurements for the Crab pulsar. This is be- 
cause the material within or around the nebula can po- 
tentially influence, or may even dominate, the observed 
scattering on some (or most) occasions. In fact, ob- 
servations have shown that the dispersion and scatter- 
ing measurements of the Crab v ary significantly with 
time (iRankin fc CounselmaiJll973t iLvne fc Thorneill975t 
llsaacman fe Rankidll977tlBackenl2000HLvne et al.ll200lh . 
The extreme cases are the so-called "scattering events" 
where the pulse-broadening time increased dramatically 
(by two orders of magnitude) and subseq uently decreased 
to no rmal over a period of several months (jLvne k Thome] 
Il975[ ). Such anomalies were interpreted in terms of an 
increase in turbulence or density within the nebular re- 
gion. Barring such exceptional cases, the measurements 
of pulse-broadening time at any given frequency still show 
large variations from one epoch to the other. For example, 
ISallmen et ail (|l999h reported the broadening time chang- 
ing from 0.28 to 1.3 ms from their Green Bank ob servations 
at 300 MHz. More recently, iPopov et all (|2006l ) measured 
a value of 43 ± 5 /is at 600 MHz, which is a fac tor two 
smaller than that reported bv lSallmen et al.l (|l999l ) at this 
frequency. Fig. [5] shows a summary of these measurements 
along with those from our observations. Incidentally, our 
measurements are ~3-5 times smaller than that predicted 
based on the published measurements and frequency scal- 
ing. This may suggest that the nebular contribution to 
the total scattering is now smaller compared to those in 
earlier observations. In other words, the bulk of our ob- 
served scattering is likely resulting from the distributed 
ISM along the line-of-sight. The best-fit frequency scal- 
ing of Td oc t/- 3 - 5 ±°- 2 derived b y us (Fig. E]) is in agree- 
ment with that of IPopov et alJ ( 2006h but is a significant 
departure from a Kolmogorov scaling (i/~ 4,4 ) favored by 
some of the earlier o bservations (jlsaacman fc Rankinlll977t 
ISallmen et al.lfl999T ). when the bulk of the observed scat- 
tering was attributed to material either inside or around 
the synchrotron nebula. 




60 B0 aoo 
Frequency (MHz) 



400 600 000 



Fig. 5. — Low frequency measurements of pulse broadening time 
from recent observatio ns of the Crab pulsa r. Observations at 23, 111 
and 600 MHz are from IPopov et al.l [(20061 ). while those at 3 MH z 
and the higher value at 600 MHz are from ISallmen et al. (1999). 
The solid line shows our best-fit scaling (r^ oc v~ 3-5±0.2^ £ Q ^. ne 
data (excluding higher values at 300 and 600 MHz), which is sig- 
nificantly shallower than the v~ 4,4 scaling expected from a Kol- 
mogorov model. Our measurements (at 200 MHz) are ~3— 5 times 
smaller than that predicted based on the published measurements 
and frequency scaling. 



In summary, our measurement of a low scattering and an 
SM estimate that is somewhat lower than that predicted 
by the NE2001 electron density model (which largely ac- 
counts for the scattering due to the distributed ISM) , sug- 
gest that the contribution from the nebula itself is proba- 
bly at its lowest. This is further supported by a reasonably 
good fit obtained with PBF2 which corresponds to a dis- 
tributed scattering geometry. If the scattering due to the 
nebula was dominant, one would have expected the mea- 
sured pulse, and the implied PBF, to have quite different 
shape given the filamentary structures of the nebular ma- 
terial. Such filaments would mean scattering structures 
with large axial ratios, which may produce PBFs that are 
of the form oc t -1 / 2 exp(— t) and possibly tr uncated beyond 
certain times fe.g. ICordes fc Laziol l200lT ) . Interestingly, 
this measurement of low scattering is also consistent with 
our recent observations at a frequency of 1400 MHz made 
with the ATC A, where we estim ate the pulse broadening 
to be 5 1 /xs (|Bhat et al.l [2007). Finally, we believe such 
a low scattering is what enabled a successful detection of 
the giant pulses with our ED system. For instance, if the 
scattering was to be dominated by the nebula, increasing 
the net pulse broadening by a factor 5 to 10, then the 
peak amplitude of the pulse will be reduced by a similar 
factor, rendering a positive detection rather difficult with 
our simple equipment. All in all, the measurement of a 
low scattering and our detection prospects are coupled, 
and perhaps it is a coincidence that our observations were 
made at an epoch when the scattering happened to be 
unusually low. 

4. GIANT PULSE DETECTION PROSPECTS WITH THE 
FULL MWA-LFD SYSTEM 

Our successful detection of the Crab giant pulses with 
the ED system opens up some exciting prospects for the 
studies of giant pulses and other fast radio transients with 
the full MWA-LFD system when it becomes operational. 
The full system will have a sensitivity that is over two 
orders of magnitude larger than the ED system, and in 
combination with a much larger bandwidth feasible for 
data recording and the instrument's unique capability to 
observe at more than one frequency simultaneously, the 
MWA-LFD will offer some brand new avenues of investi- 
gation in this arena. 

Fig. [H] shows giant pulse detection prospects of the 
MWA-LFD system, where the sensitivity curves for the 
full system are shown (middle and right panels) along 
with those for the prototype ED system (left panel). The 
sensitivity curves are dramatically different for the two 
cases, which is due to an interplay of different factors in- 
volved in the calculations. For the ED system, the sys- 
tem noise dominates over all other contributions and a 
maximum in sensitivity is expected near 200 MHz. For 
"Crab-like" pulsars that reside in "Crab-nebula-like" en- 
vironment, the sensitivity calculations take into consider- 
ation factors such as, contribution from the nebular emis- 
sion to the system noise, as well as frequency dependencies 
of the basic system performance parameters (i.e., T sys and 
G) and the sky background (T s ky). The combination of 
a decreasing T s k y and a reduced nebular contribution as 
the nebula becomes resolved at higher frequencies make 
up for the decrease in tile gain and a larger T sys expected 
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Fig. 6. — Giant pulse detection sensitivity of the MWA-LFD system. The sensitivity is expressed as maximum distance Dmax to which a 
giant pulse is potentially detectable. For the prototype system (left panel), the achievable sensitivity is dominated by the system noise itself, 
and an optimal performance is expected near 200 MHz. The middle and right panels show projected sensitivities for the full system. The 
nebular emission is important for "Crab-like" pulsars that reside in "Crab-nebula-like" environment (middle panel). The nebular contribution 
to the system noise decreases significantly above 190 MHz as the nebula is resolved by the array beam (assuming an array diameter of sal km), 
thus significantly increasing the system sensitivity at these frequencies. For field pulsars however (right panel), Dmax is largely determined 
by the sensitivity of the array itself, with the best performance expected near 200 MHz. 



at higher frequencies. This effectively yields higher sensi- 
tivity at frequencies k 200 MHz. Naturally, much higher 
sensitivities can be achieved for field pulsars with no neb- 
ular background, where D max is largely determined by the 
system performance parameters such the array gain and 
noise temperature. 

To elaborate further on our sensitivity calculations, the 
radiometer equation can be re- written for the case of gian t 
pulse detection as ( see e.g. iMcLaughlin fc Cordesl (120031 ): 
ICordes et all (120041) 1. 




where (S/N)det is the signal-to-noise achievable for a giant 
pulse of peak amplitude Sgp using a system of sensitivity 
S sys , with Av being the bandwidth over which the data are 
recorded and W p is the effective pulse width, and N po iis 
the number of polarizations that are summed. Conversely, 
for a given system with sensitivity S sys , the maximum dis- 
tance to which a giant pulse of amplitude Sgp is potentially 
detectable is given by 

Dmax= (£) '0' (Np ° iAi/Wp)i (6) 

For the MWA-LFD system, S sys is a strong function of 
the observing frequency, as both the syste m temperature 
and g ain vary significantly with frequency (jBowman et al.l 
120061) . Further, the full array will comprise 512 tiles 
spread over a region of — 1 km. Thus, both the ef- 
fective gain and resolution of the array (i.e., half-power 
width of the phased array beam pattern) are such that 
S sys is strongly influenced by the nebula. For an array 
layout extending out to 1 km, the nebula remains un- 
resolved at the lower half of the MWA-LFD 's operating 
frequency range, but becomes resolved near 200 MHz, 
where the effective resolution (beam width) becomes com- 
parable to the Crab's characteristic diameter (5. '5). The 
net system sensitivity (in units of Jy) is thus given by 
Ssys {v) — S sy so {v) + i v Scn {v), where S 

sysO IS the sys- 
tem sensitivity without the nebular contribution to the 
system noise and f„ is the beam dilution factor, given by 



f„ = Q,a{v) I&cn where fiAand ficN are the solid angle 
areas of the phased array beam and the nebula, respec- 
tively. The system sensitivity when pointed away from 
the nebula is S sys o (f) = T sys o (v)/G(u), where T sys [v) = 
T s k y (^)+T rccv (u) is the net system temperature. Further- 
more, the nebular emission is also frequency dependent 
(see § I2.3.3j) , and the integrated flux density is expected 
to change as much as by 30% over the MWA-LFD's fre- 
quency range. 

Thus, in the lower half of the MWA-LFD frequency 
range (100-200 MHz), sensitivity is largely determined by 
the combination of a decrease in T s k y and an increase in 
the tile gain, and then sharply increases near mid-way as 
the nebula becomes resolved (Fig. [5]). A giant pulse as 
bright as the brightest one in our data is thus potentially 
detectable with the full system out to a distance of — 17 
kpc at 100 MHz and out to -30 kpc at 300 MHz. Indeed 
these limits will be even larger in the case of giant-pulse 
emitters that do not reside in a nebula. In short, bright 
giant pulses are potentially detectable over a wide range 
of Galactic distances within the full range of MWA-LFD 
frequencies. 

The brightest giant pulse detected in our data (Sgp = 
50 kJy) has an energy content of — 1.85 x 10 5 kJy fis. 
Although our ED system is not capable of detecting the 
integrated emission from the pulsar, assuming a pulsar flux 
density of — 650 mJy at 400 MHz and a spectral index of 
—3.1 (e.g. iManchester et al.ll2005h . it turns out such a gi- 
ant pulse will be at least — 100 times more energetic than 
the average emission (if we assume the spectrum does not 
turn over until —200 MHz). Given this, and an expected 
two orders of magnitude improvement in sensitivity for the 
full array, it is quite likely that almost all pulses emitted 
through the giant-pulse emission mechanism (say, energy 
10 x the average emission) will be potentially detectable 
with the full system at 200 MHz. This will enable a de- 
tailed characterisation of the giant pulses and their energy 
distributions at low frequencies, a poorly studied aspect of 
the Crab giant pulses. 

It is worthwhile to examine the detection prospects for 
other known giant-pulse-emitters such as PSRs B0540- 
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69 and B1937+21. PSR B0540-69 is a young, Crab-like 
pulsar located inside the Large Magellanic Cloud (D ~ 
50 kpc), for which the giant-pulse energy distribution is 
known to follow the form P(E > E a ) = KE~ a ; where 
P is the probability of emitting a giant pulse with energy 
greater than E a (in units of t he average pulse energy), and 
X=0.26 and a ~ 1.8 from Uohnston et alj (|2004l ). The 
pulsar period is 50.35 ms and thus the strongest pulse 
emitted in 1 hr will have an energy E a ~ 235(E) where (E) 
is the average emission. Assuming giant pulses dominate 
the pulsar emission and a spectral index of —3.6 (and no 
spectral turnover down to ^200 MHz), we estimate the 
pulsar flux to be 27 mJy at 200 MHz, thus a giant pulse 
with energy 235(E) would imply 3.3 x 10 2 kJy fj,s. This 
is 56 times weaker than our brightest giant pulse, or 20 
times below the detection threshold of the ED system. 
However, with the full array (~170 times improvement in 
sensitivity), such a giant pulse should be detectable with 
S/N ~ 85. And given the KE~ a form for the energy 
distribution, we estimate ~ 50 pulses to be potentially 
detectable in 1 hr of observation. 

For the millisecond pulsar B1937+21, the giant-pulse en- 
ergies follow a simil ar power-law form wi th K=0.032 and 
a ~ 1.8 at 430 MHz (jCognard et al.lll996f) . Following sim- 
ilar arguments as above, we estimate the brightest giant 
pulse from this pulsar in 1 hr will be ~ 500 times as ener- 
getic as its average emission. While this is 8 times below 
the detection threshold of our ED system, it will be easily 
detectable (S/N ~ 212) with the full array and assuming 
a power-law index of 1.8, we estimate there to be as many 
as ~ 250 giant pulses with S/N > 10 in 1 hr duration. 
While these are only some rough estimates, and do not 
take into consideration the impact of increased scattering 
at low frequencies, it is evident that the full system offers 
exciting prospects for giant pulse studies in general. 

While the above calculations are specifically performed 
for giant pulses, the same treatment also applies to fast ra- 
dio transients of similar characteristics (i.e., time scales of 
the order of 100 [is or longer). The superb RFI-quiet envi- 
ronment of the proposed MWA-LFD site can be exploited 
to take the best advantage of this potential. Moreover, the 
instrument's capabilities such as multibeaming will serve 
as powerful discriminators of real versus spurious signals. 
Thus the full system is indeed a promising instrument for 
science related to giant pulses and similar fast radio tran- 
sients at low frequencies. 

5. CONCLUSIONS AND FUTURE WORK 

Using a simple equipment operating in the Western Aus- 
tralian outback, we have successfully detected giant pulses 
from the Crab Nebula pulsar at a frequency of 200 MHz. 
Despite a large sky background severely limiting the sensi- 
tivity achievable at such low frequencies, our system, com- 
prising just three tiles, each consisting of a 4 x 4 array of 
dipoles, detected a total of 31 giant pulses over a dura- 
tion of 3.5 hours. The measured pulse shape is signifi- 
cantly broadened by multipath scattering due to the ISM, 
and results in a degradation of the intrinsic pulse ampli- 
tude by nearly a factor of two. The deconvolution of the 
measured pulse using a CLEAN-based procedure yields a 
pulse-broadening time of 670±100 fxs for the case of a thin 
slab scattering model, and 400±50 fis for a model where 



the scattering material is uniformly distributed along the 
line of sight. The implied degree of scattering is the low- 
est that has been reported towards the Crab pulsar from 
observations made so far. In fact, our detection of giant 
pulses in such large numbers would not have been pos- 
sible but for such low level of scattering. Together with 
recent observations at low frequencies, our measurements 
of the pulse-broadening time support a shallow scaling in 
frequency (r<j oc j,- 3 - 5 ± - 2 ). With the sensitivity of our 
equipment (a gain of ^0.01 K Jy _1 and a system tempera- 
ture of ^400 K towards the Crab), pulses that are brighter 
than ^9 kJy in amplitude are easily detected. The bright- 
est giant pulse detected in our data has a peak amplitude 
of ^50 kJy and a width of «300 [is, and the equivalent 
brightness temperature is 10 31 ' 6 K, assuming a pulsar dis- 
tance of 2 kpc. 

The success of our exploratory observations underscores 
the tremendous potential the full MWA-LFD system will 
offer for studies of pulsar science in general, and for gi- 
ant pulses and fast transients in particular, at low radio 
frequencies. With a sensitivity that is over two orders of 
magnitude larger than that of the prototype equipment, 
the full system will have the capability to detect giant 
pulses over a wide range of Galactic distances. For in- 
stance, giant pulses as bright as ~50 kJy will be poten- 
tially detectable out to a distance of ^30 kpc in the case 
of Crab-like pulsars, and even further in the case of objects 
that do not reside in a nebula. In addition to enabling an 
in-depth study of Crab giant pulses at frequencies £300 
MHz, the MWA-LFD will be a promising instrument for 
a wide variety of pulsar and transient science at low fre- 
quencies. 
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